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Abstract 

A closed set of exact equations describing statistical theory of tur- 
bulent self-diffusion by multivariate-normal turbulent velocity field is de- 
rived. In doing so, we first suggest exact formulas for correlations (/i(p)/, (p')7?[f]), 
{g(j))R[f\) and {g{p)fj(p')R\i]) when the functional R[f\ is functional of 
functions /i's having multivariate-normal distribution, g and /^'s have 
joint normal distribution and zero mean values. 

Originative works of Taylor [1] and Batchelor [2] on turbulent self-difFusion 
of fluid particles in Lagrangian and Eulerian frameworks, respectively, have 
placed statistics of Lagrangian velocity field, first put forward by Taylor [1] , as 
fundamental properties in the field of statistical theory of diffusion. Despite 
many persistent efforts (for review, see e.g., Refs. [3, 4]), the theory of diffu- 
sion remains incomplete mainly due to involved closure problems [4, 5, 6] and 
lack of accurate prediction of Lagrangian statistics [7, 8, 9] in general turbu- 
lent flows. In this letter, we base our Eulerian analysis on the equation for 
a passive scalar field ilj{x,t) whose ensemble average (denoted by ()) denotes 
Green's function G for the evolution in space-time (x — t) of an arbitrary ini- 
tial probability distribution for the position of marked fluid particles in space 
[2, 5, 8]. We solve the closure problem arising in the equation for (^/'(x, t)) re- 
sulting in an appearance of Lagrangian velocity correlations. We then obtain 
equations relating Lagrangian and Eulerian correlations for the velocity field. 
These equations are exact when the fiuctuations u[ in turbulent velocity field 
Ui{x.,t) = Ui{x,t) + u^(x, i) over the mean velocity J7i(x, obey multivariate- 
normal distribution. In doing so, we first suggest three new functional formulas 
which along with the Furutsu-Donsker-Novikov functional formula [10] are used 
in our Eulerian analysis. These four functional formulas are now discussed 
briefly before presenting the Eulerian analysis. 
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Functional Formulas: Furutsu-Donsker-Novikov formula [10] suggests tliat 
for any random function fi{p) with normal or Gaussian distribution and {fi{p)) = 
0, its correlation with functional -R[f] of function /,'s can be exactly given by 

iMp)m]) = JiMp)Mp')){j^)dp', (1) 

where the integral extends over the region of arguments p in which the function 
fi is defined and ^ Sf^(^y)dp' ) represents the fimctional derivative of R with 
respect to fk- Here, we suggest another exact functional formula for higher 
order correlations {fi{p)fj{p')R[f]), written as 

{/.(i>)/,(p')fl[fi> = {mw)){m) + \l (/<(p)A(»)> <<• 

Yet another formulas which can be obtained from Eqs. (1) and (2), respectively, 



are 



^^^)m)-jmh¥)){j^)dp'. (3) 

{9{p)f,{p')R[i]) = {g{p)fj{p')}{R[t]} + \j {g{p)h{s)) (^f,{p')^^^ ds 

(4) 

where g and /j's have joint normal distribution and {g) = 0. It should be 
noted that the terms with functional derivatives in Eqs. (2) and (4) can be 
further simplified by using Eqs. (1) and (3). From now onward, these functional 
formulas given by Eqs. (1), (2), (3) and (4) are referred to as FI, F2, F3 and 
FA, respectively. 

Now, wc bricfiy describe equations which can be used for verifications of the 
functional formulas. Using Novikov's notation, expansion of R[f] around f = 
can be written as 

ii[f] = i?[0] + 1 f-- I R'^liJsu;Sn)fM---fiAsn)ds,---dSn (5) 

n=l 

where 

Jt,,...iJSl,-,Sn; <5/i,(si)dSi---5/i„(s„)rfs„ f=0 ^""^ 
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and summation over the repeated indices ii,i2- ■ - in is taken. The functional 
derivative of Eq. (5) yields 

6R[f] SR[f] 



Sfkis)ds 6fk{s)ds 

For mull 
and fi's 



1 /' /' 

^^P+X) (n- 1)! J " J ^fei2-i„(*'«2,---,Sn)/i2(s2)---/i„(s„)(is2---(is„. 

_ , (7) 

For multivariate- normal distribution for /^'s and joint normal distribution for g 



{A{s)fM---fi^{sm))=0, (8) 

n 

(A(s)/ii (Sl) ■■■fin (Sn)) = X] (^(*)/'» M){fii (Sl) • • • /ia-i (Sa-l)/i„+i (^a+l) •••/»„ 

a=l 

(9) 

where A(.s) is equal to any fi{s) or g{s) and m is an even integer. Using Eqs. 
(5)-(9), formulas F2, F3 and F4 can be verified. 
Eulerian Analysis: We consider scalar field 

V'(x,t)=<5(x-X(t)) (10) 

for the Eulerian description of the self- diffusion of fluid particle whose position 
X(t) (also denoted by components Xi{t)) at any time t is governed by the 
Lagrangian equation 

^=u.(x = X(t),i) (11) 

with initial condition X(io) = xq at time to. The governing equation for tp can 
be written as (see e.g. [5]) 

d d 

— V'(x,t) + —[i^{x,t)ui{x,t)] =0,Wt>to (12) 

with initial condition ijj{'x, to) = <5(x— xq). The ensemble average of ip multiplied 
by dx, i.e. {tp)dx = G(x, t|xo, to)<^x, with the Green's function G, denotes the 
probability that a fluid particle at (xq, to) should lie within a volume dx at later 
time t > to [5]. The ensemble average of Eq. (12) 

^(^(x, t)) + A [(^)C/,(x, t) + (Mix, t))] = (13) 

with initial condition {tp{x,to)) = S{x — Xq) poses closure problem due to un- 
known correlation {il){x,t)u'-{x,t)). In our analysis, we would need yet another 
scalar function iptiy, t — s) = (5(y — X(t — s)) with initial condition given at time t 
i.e. tpb{y,t) = (^b(y,i)> =<5(y-X(t)) =S{y-x). In fact, (■i/'6(y,t-s))dy with 
s > represents the probability that a fluid particle lies within a volume dy at 
(y, t — s) if the particle reaches (x, t). The governing equation for ipbiy, t — s), 
after introducing notation = t — s, is 

g^My^ s~) - Q^\'^b{y, s-)ui{y, s-)] = 0, (14) 
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which is obtained by using the Lagrangian equation 



—^Xi{s-) = Ui{y = X(s-), S-), yO<s<t-to (15) 
as 

describing the backward evolution of trajectory of fluid particle with known 
initial position X(t) at initial time s~ = t or s = 0. The ensemble average of Eq. 
(14) poses closure problem due to the unknown correlation {ipb{y, s~)u'i{y, s~)). 

Now we obtain exact closed expressions for the unknown correlations by 
using Fl and following a method similar to that provided earlier in the context 
of kinetic approach for two-phase turbulent flows [11]. Considering ^p[u'] as a 
functional of u-'s having multivariate-normal distribution and applying Fl, we 
obtain 

(V(x,i)uUx,t)) = - A[A,,(V,)] +7,(V>) (16) 



dxk 



where tensors 



Xki = / {4{x,tyj{t\t2))Gjk{t2\t)dt2, (17) 

Jto 

li = J' {^^^^u'jit\h))Gjr{t2\t)dt2, (18) 

and shorthand notation u'^{t\t2) is used to represent u'^ at (X(t2),i2) along the 
fluid particle path which reaches x at t. Also Gjk is governed by 

j^Gjk{t2\t) - ^E^GMt) = SjkS{t - t2) (19) 

and Gjk = Sjk for homogeneous turbulence with uniform mean velocity Uk- The 
Xki represents eddy diffusivity and ji — represents drift velocity when Eq. 
(16) is substituted in Eq. (13). We should mention that for the homogeneous 
turbulence with constant Uk, the present solution for {'ipu'^) becomes identical 
to the solution derived by Reeks [12]. Also, the present solution is valid for both 
incompressible and compressible velocity fields. Similarly, considering i/jf, as a 
functional of it^(y, s^) and applying Fl, we obtain 

{My,s-K{y,s-)) = - A[A^^(^,)] +^6^V6) (20) 



where the tensors 



AL=/ {K{y,s-yj{s-\s2))G%{s2\s-)ds2, (21) 
Jo 



l\ ^''%P ^>~\'2))Gl{s2\s-)ds2, (22) 
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and u'j{s ) represents u'^ at time = t — S2 along the particle trajectory 
which passes through y at time s~ =t — s. Also, G^fe(s2^|s~) satisfies 

-jG%{s^\s-) - ^Eh^^Glis^ls-) = Sj,S{s- - s^). (23) 

The right hand side (rhs) of Eqs. (17), (18), (21) and (22), contain unknown 
Lagrangian velocity field correlations {u^{x, t)u'j{t\t2)) , {u'j{t\t2)du[{yi,t) / dxr) , 
{u^{y,s~)u'j{s^\s2 )) and {uj{s~\s2)du^{y,s~)/dyr}, respectively. The exact 
equations relating these correlations to the Eulerian velocity field correlations 
can be derived by using the functional formulas. Using ip and ipb, the correlations 
can be written in the forms 

{0'[a']Uj{t\t2)) = J (a{„,]^(y2,s^)V6(y2,S2 )>c^y2, (24) 

(C|^,]^(S"|S^)) = j (c(^,]^(x2,t2)V'(x2,i2))rfX2, (25) 

where a|^,] = u^(x, t) or a|^,] = du[{x,t)/dxr and c|^,] = ^^(y, .s~) or c|^,] = 
9u-(y, s~)/dyr- Here in Eqs. (24) and (25), we have introduce another notations 
[a'] and [c'] for subscripts and which are equivalent to the subscripts of selected 
variables for a'^,j and Cj^,j. For example, [a'\ = i and [a'] = ir when a'^,j = Ui 

and a|^,j = dui/dxr, respectively. Also, ^Jb{'y2iS2) = 5{y2 ~ ^(■*2^)) ^-nd t2 = 
t — S2 = S2 with initial condition given at time i or s = and initial condition 
for ^(x, to) = (5(x — y) is given at time to = t — s. It should be noted that 
coupling between the equations for (ijj) and {ipb) arises due to Eqs. (24) and 
(25). We now obtain expressions for third order correlations appearing on the 
rhs of Eqs. (24) and (25). 

We introduce a few shorthand notations as (a) = (x, <), (b) = (y, s~), 
(n) = {yn,s~) and dn = dy-ndsn and t„ = t - s„ = s~ for n = 1,2, 
Now consider general form {a'^^,^u'j{y, s~)tpb{y, s~)) with tpb as a functional of 
Uj{y, s~) and by applying Fl, F2 or F3, FA depending on whether a|^,j is equal 
to u'^ or not, we can obtain expression for the third order correlation as 

(a|„,](a)^(b)V6(b)) = (a|„,p;.(b))(V6(b)) + A[a']j, (26) 

with 

92 d d 

where the tensor 

^[a']jpi=J dsi J ds3(a|„,]4(s-|sr))(w;(b)u;(s-|s3))G^p(s3 |s-)G^;(sr|s-) 

(28) 
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and {s~ \s~) is used to represent value of velocity fluctuation at (X.(t—Sn),t—Sn) 
along the fluid particle trajectory which passes through y at time t — s. It should 

be noted that < s„ < s. The expressions for remaining tensors Cl^a']jp^ ^[a']ji 
and ^[a']j in Eq. (27) arc equivalent to the rhs of Eq. (28) but with changes 

u'jih) ^^^^jp, u'jO^) and u^(b) -> respectively. It should be 

noted that substituting only the first term on the rhs of Eq. (26) into Eq. (24) 
would yield results obtained by using Corrsin's hypothesis [7]. The remaining 
non-zero terms A[a']j and B[a']j in Eq. (26) account for the correlation between 
"[a']'^ and ipb - {ipb)- 

Now, the tensors present in Eq. (27) contain Lagrangian correlations of the 
form (cj^/]Wj (s~|s~)) and for which expressions can be written similar to that of 
Eq. (25). To obtain closed expression for the third order correlation appearing 
in Eq. (25), we consider general form (cj^,jM^(a)f/'(a)) and by applying F1,F2 
or F3, FA depending on whether c|^,j is equal to u- or not, we obtain 

(c|,,p;.(a)V(a)) = (c|,,]n;.(a))(V(a)) + C[,,],. (29) 

with 

92 d d 

where the tensor 

A[c']j>i = f dh f (it3(c[e']4(*|ii))K(a)<(i|t3))Ggp(t3|i)G«(ii|i) (31) 

The expressions for remaining tensors cij[c']j"p) '"'[c'jj; and 7[c']j in Eq. (30) are 
equivalent to the rhs of Eq. (31) but with changes Uj(a) '~d^^ ^'j^^ ~^ 

and Uj (a) 'ox'dx] ' respectively. Further, the tensors present in Eq. 
(30) contain Lagrangian correlations of the form (a'^,jU^(t|t„)), similar to Eq. 
(24) and for which expressions are derived above. 

Concluding Remarks: In this letter, it has been shown that a closed set of 
exact statistical equations describing the turbulent self-diffusion by multivariate- 
normal velocity field can be obtained after solving involved closure problems 
by using Furutsu-Donsker-Novikov formula along with other three functional 
formulas. The only unknowns in the set are various Eulerian correlations of the 
velocity field. We have seen that the statistical theory of the diffusion requires 
two scalar fields ip and tpb describing the forward and backward calculation of 
the diffusion with initial conditions given at time to < t and t, respectively, for 
the evolutions of ip and ipb- Further it has been shown that the Lagrangian 
correlations appearing in the theory can be related to Eulerian correlations 
through the use of tp and tpb and solving the involved closure problems in an exact 
manner by using the functional formulas. The Eulerian analysis presented in 
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this letter is valid for both homogeneous and inhomogeneous turbulent velocity 
field. Lastly, we should mention that the functional formulas can also be used 
to derive expressions for Lagrangian correlations of fluid velocity field along 
the inertial particles which occur in the kinetic approach for dispersed phase of 
particles/droplets in turbulent flows ( e.g., see Refs. [13] and references cited 
therein) . 
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